Understanding magmatic processes and seismo-volcano source localization with multicomponent seismic arrays

Dans cette thèse, nous étudions le problème de la localisation de sources sismo-volcanique, à partir des données enregistrées par des réseaux de capteurs composés de nouveaux sismomètres à trois composantes (3C). Nous nous concentrerons sur le volcan Ubinas, l'un des plus actifs au Pérou. Nous développons une nouvelle approche (MUSIC-3C) basée sur la méthode MUSIC permetant de retourner les 3 paramètres utiles (lenteur, azimut et incidence). Pour valider notre méthodologie, nous analysons des sources synthétiques propagées en tenant compte de la topographie du volcan Ubinas. Dans cette expérience, les données synthétiques ont été générées pour plusieurs sources situées à différentes profondeurs sous le cratère Ubinas. Nous utilisons l'algorithme MUSIC-3C pour les relocaliser. Nous traitons également des données réelles provenant d'une expérience de terrain menée sur le volcan Ubinas (Pérou) en 2009 par les équipes de recherche de l'IRD-France (Institut de Recherche pour le Déveleppment), UCD l'Irlande (projet VOLUME) et l'Institut de Géophysique du Pérou (IGP). Nous utilisons l'algorithme MUSIC-3C pour localiser les événements explosifs (type vulcanien), ce qui nous permet d'identifier et d'analyser les processus physiques de ces événements, à la suite de cette analyse, nous avons trouvé deux sources pour chaque explosion situées à 300 m et 1100 m en dessous du fond du cratère actif. Basé sur les mécanismes éruptifs proposés pour d'autres volcans du même type, nous interprétons la position de ces sources ainsi que les limites du conduit éruptif impliqué dans le processus de fragmentation.In this thesis, we study the seismo-volcanic source localization using data recorded by new sensor arrays composed of three-component (3C) seismometers deployed on Ubinas stratovolcano (Peru). We develop a new framework (MUSIC-3C) of source localization method based on the well-known MUSIC algorithm. To investigate the performance of the MUSIC-3C method, we use synthetic datasets designed from eight broadband isotropic seismic sources located beneath the crater floor at different depths. The fundamental scheme of the MUSIC-3C method exploits the fact of the cross-spectral matrix of 3C array data, corresponding to the first seismic signal arrivals, provides of useful vector components (slowness, back-azimuth and incidence angle) from the seismic source. Application of the MUSIC-3C method on synthetic datasets shows the recovery of source positions. Real data used in this study was collected during seismic measurements with two seismic antennas deployed at Ubinas volcano in 2009, whose experiment conduced by volcanic teams of IRD-France (l'Institute de Recherche pour le Déveleppment), Geophysics group University College Dublin Ireland and Geophysical Institute of Peru (IGP). We apply the MUSIC-3C algorithm to investigate wave fields associated with the magmatic activity of Ubinas volcano. These analysis evidence a complex mechanism of vulcanian eruptions in which their seismic sources are found at two separated sources located at depths of 300 m and 1100 m beneath the crater floor. This implies the reproduction of similar mechanisms into the conduit. Based on the eruptive mechanisms proposed for other volcanoes of the same type, we interpret the position of this sources as the limits of the conduit portion that was involved in the fragmentation process.SAVOIE-SCD - Bib.électronique (730659901) / SudocGRENOBLE1/INP-Bib.électronique (384210012) / SudocGRENOBLE2/3-Bib.électronique (384219901) / SudocSudocFranceF

Permanent tremor of Masaya Volcano, Nicaragua' Wave field analysis and source location

International audienceThe Masaya Volcano, Nicaragua, is a basaltic caldera in a subduction zone. The permanent source of the volcanic tremor was located inside Santiago crater, at the lava lake's position and 400 m below the NE rim, and therefore corresponds to superficial magma activity. We used two tripartite arrays (90 m side), one semicircular array (r=120 m) in 1992, and two semicircular arrays (r=60 m) and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993. We used both a cross-spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field. Using the delays therefore by intersecting the back azimuth wave directions from the arrays, we could pinpoint the source. Additionally, the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor. The tremor's wave field is composed of comparable quantities of dispersed Rayleigh and Love surface waves, whose phase velocities lie in the ranges 730-1240 m/s at 2 Hz and 330-550 m/s at 6 Hz. The dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers. The resulting velocity models are similar for the northern and southern parts of the volcano. After geometrical spreading corrections, Q2Hz = 14 and Q3Hz =31 were determined along the northern linear array. The typical low velocities and low Q corresponding to the cone structure and are similar to those of other basaltic volcanoes like Puu Oo, Hawaii, and Klyuchevskoy, Kamchatka

Permanent tremor of Masaya Volcano, Nicaragua' Wave field analysis and source location

International audienceThe Masaya Volcano, Nicaragua, is a basaltic caldera in a subduction zone. The permanent source of the volcanic tremor was located inside Santiago crater, at the lava lake's position and 400 m below the NE rim, and therefore corresponds to superficial magma activity. We used two tripartite arrays (90 m side), one semicircular array (r=120 m) in 1992, and two semicircular arrays (r=60 m) and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993. We used both a cross-spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field. Using the delays therefore by intersecting the back azimuth wave directions from the arrays, we could pinpoint the source. Additionally, the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor. The tremor's wave field is composed of comparable quantities of dispersed Rayleigh and Love surface waves, whose phase velocities lie in the ranges 730-1240 m/s at 2 Hz and 330-550 m/s at 6 Hz. The dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers. The resulting velocity models are similar for the northern and southern parts of the volcano. After geometrical spreading corrections, Q2Hz = 14 and Q3Hz =31 were determined along the northern linear array. The typical low velocities and low Q corresponding to the cone structure and are similar to those of other basaltic volcanoes like Puu Oo, Hawaii, and Klyuchevskoy, Kamchatka

Exploring Arctic Transpolar Drift During Dramatic Sea Ice Retreat

International audienc

Permanent tremor of Masaya Volcano, Nicaragua' Wave field analysis and source location

International audienceThe Masaya Volcano, Nicaragua, is a basaltic caldera in a subduction zone. The permanent source of the volcanic tremor was located inside Santiago crater, at the lava lake's position and 400 m below the NE rim, and therefore corresponds to superficial magma activity. We used two tripartite arrays (90 m side), one semicircular array (r=120 m) in 1992, and two semicircular arrays (r=60 m) and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993. We used both a cross-spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field. Using the delays therefore by intersecting the back azimuth wave directions from the arrays, we could pinpoint the source. Additionally, the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor. The tremor's wave field is composed of comparable quantities of dispersed Rayleigh and Love surface waves, whose phase velocities lie in the ranges 730-1240 m/s at 2 Hz and 330-550 m/s at 6 Hz. The dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers. The resulting velocity models are similar for the northern and southern parts of the volcano. After geometrical spreading corrections, Q2Hz = 14 and Q3Hz =31 were determined along the northern linear array. The typical low velocities and low Q corresponding to the cone structure and are similar to those of other basaltic volcanoes like Puu Oo, Hawaii, and Klyuchevskoy, Kamchatka

Permanent tremor of Masaya Volcano, Nicaragua' Wave field analysis and source location

International audienceThe Masaya Volcano, Nicaragua, is a basaltic caldera in a subduction zone. The permanent source of the volcanic tremor was located inside Santiago crater, at the lava lake's position and 400 m below the NE rim, and therefore corresponds to superficial magma activity. We used two tripartite arrays (90 m side), one semicircular array (r=120 m) in 1992, and two semicircular arrays (r=60 m) and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993. We used both a cross-spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field. Using the delays therefore by intersecting the back azimuth wave directions from the arrays, we could pinpoint the source. Additionally, the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor. The tremor's wave field is composed of comparable quantities of dispersed Rayleigh and Love surface waves, whose phase velocities lie in the ranges 730-1240 m/s at 2 Hz and 330-550 m/s at 6 Hz. The dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers. The resulting velocity models are similar for the northern and southern parts of the volcano. After geometrical spreading corrections, Q2Hz = 14 and Q3Hz =31 were determined along the northern linear array. The typical low velocities and low Q corresponding to the cone structure and are similar to those of other basaltic volcanoes like Puu Oo, Hawaii, and Klyuchevskoy, Kamchatka

Permanent tremor of Masaya Volcano, Nicaragua' Wave field analysis and source location

International audienceThe Masaya Volcano, Nicaragua, is a basaltic caldera in a subduction zone. The permanent source of the volcanic tremor was located inside Santiago crater, at the lava lake's position and 400 m below the NE rim, and therefore corresponds to superficial magma activity. We used two tripartite arrays (90 m side), one semicircular array (r=120 m) in 1992, and two semicircular arrays (r=60 m) and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993. We used both a cross-spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field. Using the delays therefore by intersecting the back azimuth wave directions from the arrays, we could pinpoint the source. Additionally, the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor. The tremor's wave field is composed of comparable quantities of dispersed Rayleigh and Love surface waves, whose phase velocities lie in the ranges 730-1240 m/s at 2 Hz and 330-550 m/s at 6 Hz. The dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers. The resulting velocity models are similar for the northern and southern parts of the volcano. After geometrical spreading corrections, Q2Hz = 14 and Q3Hz =31 were determined along the northern linear array. The typical low velocities and low Q corresponding to the cone structure and are similar to those of other basaltic volcanoes like Puu Oo, Hawaii, and Klyuchevskoy, Kamchatka

Automatic classification of seismo-volcanic signals at La Soufrie;re of Guadeloupe

International audienceSeismic activity at La Soufrière volcano of Guadeloupe is composed of various transient signals, which are classified manually by the Observatoire Volcanologique et Sismologique de Guadeloupe (OVSG-IPGP) considering waveforms recorded at several stations. Although five main types of signals are recognized in the data analysis by the observatory (Moretti et al., 2020), only three main classes readily distinguishable on seismic traces during the daily analytical protocol have been catalogued: Volcano-Tectonic events, Long-Period events and Nested events, each related to a distinct physical process. Automatic classification of seismo-volcanic signals of La Soufrière was performed by using an architecture based on supervised learning, available at github.com/malfante/AAA. Seismic waveforms are transformed into a large set of features (34 features for each representation domain) computed from three representation domain of the signal (time, frequency, quefrency). The resulting vectors of features are then used for the modeling. We are using the Random Forest Classifier algorithm from the scikit-learn library. At first, we trained the model with the dataset given by the OVSG consisting of 845 available labeled events (542 VT, 217 nested and 86 LP) recorded in the period 2013-2018. We obtained an average classification rate of 72 %. We determined that the VT class includes a variety of signals covering the LP, Nested and VT classes. Reviewing in details the waveforms and the spectral characteristics of the signals belonging to the 3 classes we then introduced Hybrid events and also defined a monochromatic class (so-called Tornillo) of LP signals, thus matching the full description of signals provided in Moretti et al. (2020). Then, using the new information, a new model was trained with 5 classes and tested. We obtained a much better classification average rate of 84 %. The classification is excellent for Nested events (93 % of accuracy and precision) and Tornillo events (93% of accuracy and precision). The classification of VT events (90% accuracy, 89% precision) and LP events (86% accuracy, 82% precision) were also very good. The most difficult class to recognize is the Hybrid class (64 % accuracy, 69 % precision). Hybrid events are often mixed with VT and LP events. This may be explained by the nature of this class and the physical process that includes both a fracturing and a resonating component with different modal frequencies. Machine learning is a powerful tool to handle large datasets. From a dataset built manually, the processing we applied allowed to obtain a reliable automatic classification by refining class definitions. This has important implications for observatory data processing during unrest and eruptive activity

La crisis sísmica en el volcán Irazú en 1991 (Costa Rica)

Irazú is a stratovolcano located in the Central Volcanic Range, 20 km NE from San José, capital of Costa Rica. This volcano has a very important record of activity in historic times. The last eruptive period was in 1963-1965. A seismic swarm was detected at lrazú, starting in January 1991, after an important earthquake (MD 5,8) recorded 50 km away from the volcano. This swarm lasted for several weeks and included events of magnitudes strong enough to be felt by the people in San José. After several weeks of repose, another important earthquake on April 22, 1991, located 80 km away from Irazú reactivated the seismic swarm near the volcano, with a rapid increase in the number of events. Three main types of seismic signals have been identified at a) volcano-tectonic earthquakes (A-type), b) low frequency earthquakes (B-type) and c) tremor. The vulcano-tectonic activity was more important during the half of the year 1991 and low frequency volcanic events were predominant after that. Volcano-tectonic earthquakes were distributed mainly l to 5 km SSE of Irazú crater, with depths frorn 5 to 10 km, B-type events concentrated within a radius of l km around and the summit. The energy budget estimated for the present seismic crisis was in the order of 1018 ergs. Spectral analysis performed for the low frequency events showed frequency picks of 1.2, 2.9 and 2.3 Hz; tremor spectral content was in the range of l to 4 Hz. The l991 crisis can be interpreted according to two different hypothesis (or combination of both): a) A magmatic movement associated to a new magma intrusion or residual magma from the last emptive period, b) Unstabilization of the volcanic system by regional tectonic earthquakes. The lack of clear evidences for a new shallow magmatic intrusion supports the hypothesis of disturbance in the volcanic tectonic ond hydrothermal system by the strong tectonic earthquakes in december 22, 1990 and April 22, 1991.El Irazú es un estratovolcán localizado en la Cordillera Volcánica Central. Este volcán tiene un importante registro de actividad en tiempos históricos. El último período eruptivo fue el de 1963-1965. Un enjambre sísmico fue detectado en el Irazú a principios de enero de 1991 , luego del terremoto de Piedras Negras del 22 de diciembre de 1990 de magnitud MD 5,8 Este enjambre se prolongó por varias semanas y tuvo eventos sísmicos de magnitudes hasta 4,3, que pudieron ser sentidos inclusive en el Valle Central. Luego de varias semanas de reposo, otro terremoto importante, el Terremoto de Limón del 22 de abril de 1991, reactivó el enjambre sísmico cerca del volcán, con un rápido incremento en el número de eventos. Tres tipos principales de señales sísmicas han sido identificadas en el Irazú: a) sismos volcano-tectónicos (tipo A), b) sismos de baja frecuencia (tipo-B) y c) trémores. La actividad volcano tectónica fue más importante durante la primera mitad del año 1991, y los sismos volcánicos de baja frecuencia fueron los que predominaron luego. Los epicentros de los sismos volcano-tectónicos se distribuyeron principalmente entre  l y 5 km al SSE del cráter del Irazú, a profundidades entre 5 y l0 km. Los eventos tipo-B se concentraron principalmente dentro de un radio de l km alrededor y bajo el cráter. La energía sísmica Iiberada durante esta crisis se ha estimado en 1018 ergios. Análisis espectrales realizados para eventos de bajo frecuencia mostraron picos de frecuencias de 1,2, 2,9 y 2,4 Hz. El espectro de frecuencias de los trémores está contenido dentro del rango de l a 4 Hz. La crisis del Irazú de 1991 puede ser interpretada de acuerdo con dos hipótesis principales, o combinación de ambas: a) un movimiento magmático asociado a una nueva intrusión magmática o magma residual del último período eruptivo. b) desestabilización del sistema volcánico por sismos tectónicos regionales. La falta de evidencias claras para una nueva intrusión magmática somera favorecen la hipótesis de una desestabilización del sistema tectónico e hidrotermal del volcán a causa de los terremotos de Piedras Negras y Limón-Telire

Low-frequency bursts of horizontally polarized waves in the Arctic sea-ice cover

International audienceWe report the detection of bursts of low-frequency waves, typically f = 0.025 Hz, on horizontal channels of broadband seismometers deployed on the Arctic sea-ice cover during the DAMOCLES (Developing Arctic Modeling and Observing Capabilities for Long-term Environmental Studies) experiment in spring 2007. These bursts have amplitudes well above the ambient ice swell and a lower frequency content. Their typical duration is of the order of minutes. They occur at irregular times, with periods of relative quietness alternating with periods of strong activity. A significant correlation between the rate of burst occurrences and the ice-cover deformation at the ~400 km scale centered on the seismic network suggests that these bursts are caused by remote, episodic deformation involving shearing across regional-scale leads. This observation opens the possibility of complementing satellite measurements of ice-cover deformation, by providing a much more precise temporal sampling, hence a better characterization of the processes involved during these deformation events